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Abstract 

Broadband implementations of time-optimal geodesic pulse elements are intro- 
duced for the efficient creation of effective trilinear coupling terms for spin sys- 
tems consisting of three weakly coupled spins 1/2. Based on these pulse elements, 
the time-optimal implementation of indirect SWAP operations is demonstrated ex- 
perimentally. The duration of indirect SWAP gates based on broadband geodesic 
sequence is reduced by 42.3% compared to conventional approaches. 



1 



1 Introduction 



In the absence of relaxation, nuclear magnetic resonance (NMR) experiments consist of a 
sequence of unitary transformations of the density operator representing the spin system 
of interest (QJ). An important goal of both theoretical and practical interest is the design 
of pulse sequences that can generate a desired unitary transformation as fast as possible, 
in order to reduce losses due to relaxation. This poses the problem of time-optimal control 
of quantum systems (013111), which is of interest for coherent spectroscopy in general, as 
well as for quantum information processing 

Here, we focus on the realization of effective propagators that correspond to the action 
of an effective Hamiltonian with trilinear coupling terms which simulate a three-spin in- 
teraction in spin chains consisting of three weakly coupled spins 1/2 (Ising coupling). In 
the context of NMR polarization-transfer experiments, methods for creating such effec- 
tive trilininear Hamiltonians have been developed and used for many years (jHlIZI)- One 
approach to create such effective Hamiltonians is based on the decoupling of certain inter- 
actions during the pulse sequence (jSHSUHl)- An earlier, more efficient approach does not 
rely on decoupling (jHl Ej) . Recently, further improved sequences were derived which 
also avoid decoupling. Even larger time savings are possible by using geodesic pulse se- 
quences ^ that can be shown to be time-optimal. However, so far the geodesic pulse 
sequences were based on weak rf pulses which severely limited the range of frequency 
offsets for which the experiment is functional. In order to pave the way for practical 
applications of these time-optimal pulse elements, we developed broadband versions of 
the time-optimal geodesic sequence. Time-optimal indirect SWAP operations (3) were 
realized experimentally in a three-spin system, demonstrating the superior performance 
of the new sequences. 
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2 Theory 



We consider a chain of three heteronuclear spins with couphng constants J12 = t/23 = J ^ 
Ji3 = and offsets i^i, z/2, and z/3. In a multiple-rotating frame (T), the corresponding free 
evolution Hamiltonian Ti^ is in general given by 



Ho = Hc + Hoffj [1] 

with the coupling term 



K = 2TiJhJ2z + 27r [2] 

and the offset term 



Hoff = ^HUji^ + 2TXU2hz + 27rz/3J3^. [3] 

The same Hamiltonian is valid if e.g. the first and third spins are homonuclear and the 
second spin is heteronuclear {vide infra). 

Many applications in NMR spectroscopy dEllZI) and NMR quantum computing (|Tn| ITT | IT2|l 

require unitary transformations of the form 

Uaf3-y{ii) = exp{-i 2-71 K hahfih^] , [4] 

where a, /3, 7 can be x, y or z. The time r required to realize such a propagator depends 
on the pulse sequence and is a function of k (see Table 1). The propagator Uaf3'y{n) can 
also be expressed as 

Uap-^{K) = exp{-i t{k) HafB-y}, [5] 

where Ti-ais-y corresponds to an effective trilinear coupling Hamiltonian of the form 

'HalS'y = 27rJejQF(K) hahfih-y, [6] 
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and the effective trilinear coupling constant Jeff{i^) is defined by 

Jed-) = [7] 

For practical applications, t{k) should be as short as possible and hence the effective 
coupling constant Jeff{i^) and the scaling factor 

should be as large as possible. For < k < 1, the theoretical limit t*{k,) for the minimum 
time required to create a propagator Wq,^^(k) is given by © 



JkU - k) 

r*{^) = ' \j \ [9] 
which corresponds to a maximum possible scaling factor 

s*{k) = [10] 
y'/€(4 - k) 

It is sufficient to consider t*{k,) and s*{k,) only for < /t < 1, because T*{2n±K,) = t*{k,), 
where n is an arbitrary integer (jHjl. 

Schematic pulse sequences corresponding to four different approaches for the creation of 
^zzz{i^) are shown in Fig. 1. These sequences can be further streamlined by reducing the 
number of pulses using well known rules {vide infra). 



Sequence A with duration 



TAi^) = [11] 



is based on the identity (jS} 

Uzzz{t^) = exp{-i7r/i^/2^} exp{-i7r k hyhz] ^^.^{^Tihzhx]- [12] 

Equivalent sequences (jH)) with the same duration ta{k) but with less pulses can be con- 
structed based on the identity 

U^^k) = Va exp{-i7r k hzhz] V^^- [13] 



4 



with 



TT TT 

Va = exp{-i-/2^} exp{-mIiJ2z} exp{-i-/2j^}. 



Equivalant sequence with the same duration ta{k) can also be constructed using 
operations 



Sequence B with duration 



J 

is based on the identity (0 El 



TT 

l^zzzili) = Vb exp{-i - K l2x} 



with 



TT 

Vb = exp{-i - l2y} exp{-i n {h.hz + hzhz)}- 



Sequence C with duration 

1 + K 

is based on the identity Q 

TT 

i^zzziii) = Vc exp{-i TT K {IiJ2y + hyhz)} V^^ exp{i - K hz}. 



with 



TT 

Vc = exp{-i - {I1J2X + hxhz)}- 



Finally, the time-optimal sequence D with duration 



k(4 - k) 



2J 

is based on the identity (jSI) 



1^zzz{k) = VdW exp{-i7rA//t(4 - K){hJ2z + hzhz) + ivr(2 - k)I2x} 



with 



Vd = exp{-i - l2y} 



[23] 



and 



W = exp{-i7r(2 - 



[24] 



The durations t{k,) and scahng factors s(k) = Jeff{n)/J of sequences A-D are shown 
in Fig. 2 and are summarized in Table 1. In all four sequences, /2-selective pulses are 
required. In addition, Ji-selective and /3-selective 180° pulses are required in sequence A 
for selective decoupling of Ju and J23 during parts of the pulse sequence. In contrast, 
sequences B-D do not require such Ji-selective or /3-selective pulses, which simplifies the 
experimental implementation of these pulse sequences if spins Ji and I3 are homonuclear. 
These sequences are also suitable for applications, where spins Ji and Is are equivalent, as 
in I2S spin systems (jHllZI)- Furthermore, the fact that decoupling is avoided in sequences 
B-D makes these experiments considerably more efficient than sequence A {vide infra). 
Sequence B is a straight-forward generalization of a well-known pulse sequence developped 
initially (jHlIZI) for the special case of k = 1. Sequence C, which has been proposed recently 
0, is more efficient than sequence B. The geodesic pulse sequence (sequence D) has the 
shortest possible duration for all values of k ©, see Fig. 2 A (top panel). For k 0, 
the duration of the geodesic pulse sequence approaches 0, in contrast to sequences A-C. 
Fig. 2 (middle panel) shows the scaling factors s(k) and Fig. 2 (bottom panel) shows the 
relative scaling factors s{k)/ n compared to the scaling factor of sequence B (sb = k). For 
K = 1, the scaling factor s of the geodesic sequence is 73.2% larger compared to sequence 
A and about 15.5% larger compared to sequences B and C. As k approaches 0, the scaling 
factor s of the geodesic pulse sequence becomes infinitely larger than the scaling factors 
of sequences A-C. For example, for k = 0.01 the scaling factor s of the geodesic sequence 
is already about a factor of 10 larger compared to sequences A and B and about a factor 
of 5 larger compared to sequence C. 
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The basic pulse sequences shown in Fig. 1 only create the desired unitary transformations 
Uzzz{i^) if all spins are on-resonance in a multiple- rotating frame, i.e. if 'Hojj = (c.f. Eq. 
3). However, for most practical applications, a finite offset range must be covered by the 
pulse sequences. Broadband versions of sequence B can be found in the literature (jEllZI)- 
Broadband versions of sequences A (jHJ and C can be created in a straight-forward way 
by inserting additional vr pulses in the existing delays to refocus chemical shift evolution. 
For example, a broadband version of sequence C is shown in Fig. 3. The robustness of the 
broadband sequence with respect to rf inhomogeneity and offsets can be further improved 
by using the x, —x, —x, x cycle flT^ for the phases of the four vr pulses which are applied 
to spins Ji and I3. 

Although no delays exist in the ideal geodesic pulse sequence shown in Fig. ID, delays 
can be introduced by replacing the weak pulse with amplitude z/^, = (2 — k)J/ ^ k(4 — k), 
duration t*{k) = ^ k(4 — k)/ (2 J), and flip angle = u^tw'^'^ = (2 — '«)7r by n hard 
pulses with flip angle a^/n and n delays of duration A = T*{K)/n. In the limit of n — cx), 
this DANTE-type (Delays alternating with Nutations for Tailored Excitation) ()13|) pulse 
sequence creates the same effect as the weak pulse and with the same limited bandwidth. 
For the present application, the DANTE-type sequence approaches the ideal sequence if 
A ^ 1/J. By inserting tt pulses in the delays of the DANTE sequence (c.f. Fig. 4), a 
broadband version of the geodesic pulse sequence can be created. As shown in Fig. 4, 
the robustness of the broadband geodesic sequence with respect to rf inhomogeneity and 
offset can be improved by using cycles or supercycles such (HH) for the 

phases of each set of four 180° pulses. 

In addition to applications in polarization transfer experiments (jElIZI), propagators corre- 
sponding to trilinear effective coupling terms are useful in the field of quantum information 
processing. For example, so-called A2 gates ()T^ can be implemented efficiently based on 
Uzzz{i^) for K = 1. Here, we focus on the implementation of SWAP operations (|TB| IT7| ITH |l 
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that make it possible to exchange arbitrary spin states of two spins in a couphng network. 
For weakly coupled spins such as in the spin system defined in Eq. (2), a direct SWAP 
gate such as SWAP{1,2) or SWAP{2,?)) between directly coupled spins Ji and I2, or 
between I2 and /a has a minimum duration of ^ 

TSWAP(1,2) = TSWAP{2,?,) = 3/(2 J). [25] 

An indirect SWAP operation SWAP{1, 3) between spins Ii and 1^, which are not directly 
coupled, can always be realized based on the following combination of the direct SWAP 
gates SWAP{1,2) and SWAP{2,3y. 

SWAP{1, 3) = SWAP{1, 2) SWAP{2, 3) SWAP{1, 2) [26] 
with an overall duration 

^SWAP{1,3) = 2 TswAP(l,2) + T5WAP(2,3) = 9/ (2J). [27] 

However, the time-optimal realization of the indirect SWAP operation SWAP{1, 3) has a 
duration of only (jSJ 

^W(i,3) = 3r*(l) = 3v^/(2J) [28] 

and hence requires only 57.7 % of the duration t|'^^p(i_3) of the conventional sequence 
(c.f. Eq. 27). This approach is based on the time-optimal realization of propagators 
Uafs^i^) (c.f. Eq. 5), which create the desired indirect SWAP{1,3) gate by the following 
sequence of operations (jHj): 

TT 

^SWAPil,3) = l^zzz{^) Uyzyi^) l^xzx{^) exp{i-l2z). [29] 

Note that all terms in Eq. 29 mutually commute and hence in experimental implementa- 
tions the order of the corresponding pulse sequnce elements is arbitrary. Based on pulse 
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sequence elements for the realization of Uzzzii^), the sequence of propagators in Eq. (29) 
can be realized in a straight-forward way by the pulse sequence shown in Fig. 5 C. The fi- 
nal 90° 2 rotations can either be implemented by a composite pulse such as 90° 90° 90° ^,or 
by adjusting the phases of all following pulses and of the receiver (I19|l . In general, z 
rotations (by angle can be implemented by an additional phase shift (by angle — of 
all following r.f. pulses that are applied to this spin and of the receiver phase for this spin 

(Uni). 

3 Experiments 

In order to test the performance of the new geodesic pulse sequences, we used the spin 
system of the amino moiety of [^^N]-acetamide as a model system (see Fig. 6) that 
corresponds closely to the model Hamiltonian Hq = He + Hajj defined in Eqs. (1-3). 
[^^N]-acetamide (Chemotrade GmbH) was dissolved in DMS0-d6 and all measurements 
were performed on a Bruker 600 MHz DMX spectrometer (Bruker Analytik GmbH) at a 
temperature of 298 K. Here, spins Ji and Is correspond to the amino protons, whereas 
h corresponds to the ^^N spin with = 88.8 Hz ^ J23 = 87.3 Hz > J13 = 2.9 Hz. 
Additional ^^J^H,^ H) and ^J{^H,^^ A^)couplings (0.7 Hz and 1.2 Hz) of h and h to the 
methyl protons of [^^N]-acetamide are about two orders of magnitude smaller than the 
^J{^H,^^N) couplings. 

Compared to a fully heteronuclear spin system, the relatively small frequency difference 
Ai/i3 = 358 Hz of the amino protons (spins Ji and I3) makes it difficult to apply short 
selective pulses to spin Ji that do not affect spin J3 (and vice versa) as required in sequence 
A (and also for its broadband implementation using additional refocussing pulses). In our 
experiments, we implemented spin-selective proton pulses by a combination of hard pulses 
and delays. For example, if spin Ji is irradiated on resonance, a selective 180°(/i) pulse 
can be implemented by the pulse sequence element 
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90°(/i,/3) - 6- 180°(/2) - 6- 180°(/2)90°(Ji,/3), 
where 5 — 2/(4Ai/i3) = 698 //s. Similarly, a selective 180° (/a) pulse can be implemented 
by the same pulse sequence element if the last 90° (/i, I3) pulse is replaced by 901j.(/i, I3). 

Based on broadband versions of sequences A, C, and D, we implemented the sequence 
shown in Fig. 5 C, which reahzes a SWAP{1,3) operation for k, — 1. Note that the 
sequence of Fig. 5 C still allows for a variation of n, which makes it possible to test 
the theoretically expected k dependence of the sequences {vide infra). For k = 1, the 
sequences for the SWAP{1, 3) gate were successfully tested for a large number of initial 
states of the spin system. Three illustrative examples are presented in Fig. 7, where 
spectra of the amino protons (spins Ji and Is) of [^^N]-acetamide are shown. Left 
(A-C) and right (A'-C) spectra reflect the states before and after an indirect SWAP{1, 3) 
operation, repectively. The initial spin states were prepared to be p(0) = hx (c.f. Fig. 
7A), (B) p(0) = 2hJ2, (c.f. Fig. 7 B), and (C) p(0) = + ^hJ^. (c.f. Fig. 7 C). As 
expected, the states of the two proton spins (spins /i and I3) are swapped for arbitrary 
initial states. 

In order to compare the durations and k dependence of the indirect SWAP{1, 3) sequences, 
we measured the efficiency of inphase transfer from Ji^ to I^^ (c-f- Fig. 7 A): For an initial 
density operator of p(0) = Iix, we defined the transfer efficiency 7713 (r) for a given pulse 
sequence of duration r as 

where {lix){0) is the initial expectation value of Iix and {I3x){t) is the expectation value of 
Isx after the pulse sequence. The corresponding experimental values of 77 were determined 
by dividing the integral of the spin 73 multiplet in the final spectrum by the integral of 
the spin Ii multiplet in the initial spectrum. 

Fig. 8 summarizes the theoretical and experimental curves of the transfer efficiency fjisir) 
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based on broadband versions of sequences, A, C, and D. In the experiments, the parameter 
K (c.f. Eq. 4) was varied in the range < k < 2. Fig. 8 A shows the theoretical r 
dependence of the transfer efficiency r^is for the pulse sequences of Fig. 1 A, C, and D, 
assuming ideal spin-selective hard pulses without rf inhohogeneity and an isolated, ideal 
three-spin system. All spins are assumed to be on-resonance (Tioif = 0) in a multiple 
rotating frame with J12 = J23 = 88 Hz and J13 = Hz. As expected (c.f. Table 1), 
transfer efficiencies of 7713 = 1 (corresponding to a complete SWAP operation) are found 
for r = 51.1 ms (sequence A), 34.1 ms (sequence C), and 29.5 ms (sequence D). 
More realistic values of the transfer efficiency 7713 (t) to be expected for our model sys- 
tem were obtained by simulating the time evolution of the density operator during the 
broadband pulse sequences for the actual coupling network (using the experimentally de- 
termined coupling constants and frequency offsets) and taking into account experimental 
pulse sequence parameters (see Fig. 8 B). Nominal rf amplitudes of 35.7 kHz and 5.5 kHz 
were assumed for ^H and ^^N pulses, respectively (corresponding to 90° pulse durations 
of 7 fis and 45 /is). The effects of rf inhomogeneity were taken into account by assuming 
a Gaussian distribution of the rf amplitudes with a full width at half hight of 10% (jSOI- 
Relaxation effects were not included. The simulated ?7i3(t) curves in Fig. 8 B qualita- 
tively match the ideal curves shown in Fig. 8 A. In particular, the position of the maxima 
appear at very similar pulse sequence durations r. However, the amplitude of the ?7i3(t) 
curves is decreased due to the effects of experimental imperfections. 
In Fig. 8 C, experimentally determined transfer efficiencies rji^lr) are shown for the three 
pulse sequences. A reasonable match is found between experimental and simulated curves. 
The experimentally determined bandwidth covered by the broadband geodesic sequence 
was about 3.5 kHz for ^H and 2.5 kHz for ^^N for the given pulse sequence parameters. 
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4 Discussion 



Broadband versions of a new class of pulse sequences for the simulation of trilinear cou- 
pling terms were developed. Using the amino group of [^^N]-acetamide as a model system, 
the theoretically predicted properties (;3| of the new sequence C and of the time-optimal 
geodesic sequence D were verified and efficient exchange of the spin state of indirectly cou- 
pled spins was demonstrated. It is expected that the new broadband pulse sequences will 
find applications both in quantum information processing and in coherent spectroscopy. 
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Table 1: Pulse sequence durations r and scaling factors s — Jeff{K)/J of effective trilinear 
coupling constants. 
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B 


c 


D 




2+K 

2J 


1 
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1 + K 

2J 


-y/K(4-K) 
2J 


r(l) 


1.5 

J 


1 
J 


1 

J 


V3 ^ 0.866 
2J ~ J 


s{k) 


2k 

2+K 


K 


2k 

l + K 


2k 

^/K{4-K) 


s{l) 


1 0.666 
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1 


75 ^ 1-155 


'TSWAP{1,3){J) 


4.5 
J 


3 
J 


3 
J 


3a/3 ^ 2.598 
2J ~ J 



''■5WAP(i,3)(88Hz) 51.1 ms 34.1 ms 34.1 ms 29.5 ms 
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-y 



(2^-|)-x 



(2it-f)^y 



Figure 1: Schematic representation of four basic (narrowband) pulse sequences for the cre- 
ation of a propagator Uzzz{,i^) = exp{— i 9 luhzhz] with 6 = 271^. If not exphcitely spec- 
ified otherwise, narrow and wide vertical bars represent spin-selective 7r/2 and vr pulses, 
respectively. (A) Example of a pulse sequence based on selective decoupling (01313), (B) 
conventional pulse sequence without decoupling (jHllZI), (C) improved sequence without 
decoupling Q, (D) time-optimal geodesic pulse sequence with r* = J k{A — k)/2J Q. 
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Figure 2: Durations r(K) (top panel), scaling factors s = Jeff{K)/J = k/Jt{k) (middle 
panel), and relative scaling factors s/sb = s/k (bottom panel) of the four basic pulse 
sequences A-D shown in Fig. ^ 
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Figure 3: Broadband version of sequence C shown in Fig. ^ 
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Figure 4: Broadband version of the geodesic sequence shown in Fig. ID for the time- 
optimal implementation of Uzzz{i^)- The pulse sequence element in brackets has a duration 
of 4A and is repeated m times. The narrow vertical bars in the bracket correspond to 
hard pulses with flip angles /n and A = r* /n with n = Am. 
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Figure 5: Pulse sequences implementing the indirect SWAP operation Uswap{i,2,) for = 
1. (A) Schematic implementation according to Eq. 29, (B) equivalent implementation 
based on pulse sequence elements (c.f. Fig. 1) that create the propagator Uzzz{,k), (C) 
streamlined pulse sequence with a minimum number of 90° pulses. 
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Figure 6: The model system (top) consisting of a chain of three coupled spins 1/2 with 
Ji2 = J23 = J and J13 = is approximated by the spins of the amino moiety (printed in 
boldface) of [^^N]-acet amide (bottom). 
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Figure 7: spectra of the amino protons (spins Ii and I^) of [^^N]-acetamide before 
(A-C) and after (A'-C) an indirect SWAP{1, 3) operation based on Fig. 5 C and the 
broadband sequence shown in Fig. 3 for the creation of Uzzzi^) for k = 1. The initial spin 
states were prepared to be (A) p(0) = (B) p(0) = 27ij.72z, and (C) p(0) = hx+'^hzhx: 
respectively. 
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Figure 8: Transfer efficiency 7713 (r) (c.f. Eq. 30) based on broadband versions of sequences, 
A (dotted curves), C (dashed curves), and D (solid curves). (A) Theoretical curves 
assuming an ideal spin system (see Fig. 6 top) and ideal rf pulses, (B) simulations based 
on the coupling constants of [^^N]-acetamide and assuming finite pulse durations and 
realistic rf inhomogeneity (see text), (C) experimental transfer curves. 
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